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INTRODUCTION 


Assessements  by  several  in-depth  workshops  have  conclusively  affirmed  the 
importance  of  continued  research  on  3-D  Shock  Turbulent  Boundary  Layer 
Interaction  ("SBLI“)  to  the  aerodynamical  engineering  of  both  external  flow 
fields  on  hypersonic  vehicles  such  as  the  NASP  and  the  internal  flows 
within  air-breathing  engine  inlets.  This  conclusion  has  been  further 
emphasized  recently  by  a  strong  surge  of  interest  in  the  accurate 
prediction  of  inlet  flows  for  air  breathing  engines  on  very  high  Mach 
number  hypersonic  flight  vehicles. 

The  present  research  involved  a  study  of  3-D  interactions  including  flow 

separation  from  nonaidabatic  surfaces,  wherein  we  sought  to  develop  a 

generalized  interactive  "wall  sublayer"  model  of  the  laminar  plus 

turbulent  disturbance  shear  stress  physics  that  embraces  the  combined 

presence  of  streamwise  flow  reversal,  cross-flow  and  heat  transfer. 

According  to  the  testimony  of  both  experimental  and  CFD  workers,  such  a 

"super  wall  layer"  model  would  not  only  provide  valuable  insight  to 

observed  interactive  skin  friction  line  behaviour  and  its  Mach  number 

dependence,  but  would  also  greatly  enhance  computational  efficiency  as  an 

"underlay"  to  a  computer  code  (say  a  3-D  Euler  code)  that  'numerically 

solves  all  the  remaining  outer  portion  of  the  flow  field.  Indeed,  the 

need  for  providing  a  careful  treatment  of  the  near-wall  region,  because  of 

the  great  difficulty  and  expense  in  providing  suff iciently-f ine  streamwise 

as  well  as  normal  numerical  mesh  spacing  at  high  Reynolds  numbers,  has  •**-“'*»* 

been  stressed  by  a  number  of  experts.  Our  research  originally  proposed' tj  tan 

has  three  main  technical  goals  within  the  scope  of  a  three  year  effort: 

(1)  The  development  of  such  a  generalized  wall  layer  model  in  the  form  of _ 
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a  set  of  first  order  partial  differential  equations  that  globally 
characterize  the  Interactive  dlsplacenent  thickness,  skin  friction  and 
heat  transfer  due  to  the  near-wall  viscous/turbulence  shear-disturbance 
physics;  (2)  Incorporation  of  this  model  as  a  wall  layer  underlay  in  a 
suitable  3-D  Euler  (or  Navier  Stokes)  code,  followed  by  validation  with 
data  from  several  basic  experimental  studies  of  shock/boundary  layer 
interactions  on  isolated  single  fins  and  ramps;  (3)  Application  of  the 
resulting  composite  code  to  the  prediction  and  analysis  of  double  fin^ 
("crossing  shock")  interactions  typical  of  the  actual  inlet  flow  field  in 
air  breathing  engine  inlets  (Fig  1). 

Owing  to  the  abrupt  and  unexpected  termination  of  the  supporting  grant 
after  only  one  year,  only  part  of  goal  (2)  and  none  of  the  work  toward 
goal  (3)  could  be  carried  out.  This  report  describes  the  progress  that  was 
achieved  regarding  goals  (1)  and  (2). 

SUMMARY  OF  PROGRESS 

At  high  Reynolds  number,  it  has  been  experimentally  and  computationally 
we 11 -confirmed  that  the  3-D  interactive  physics  is  vertically  organized 
into  the  triple-deck  type  of  structure  illustrated  in  Fig.  2,  which 
contains  at  its  bottom  a  thin  inner  shear  disturbance  deck  which  provides 
the  appropriate  matching  boundary  condition  at  the  bottom  of  the  overlying 
middle  deck  and  contains  much  of  the  interactive  displacement  thickness 
effect  plus  all  of  the  important  skin  friction  and  heat  transfer 
disturbances  and  underlying  flow  separation  zones.  It  is  the  detailed 
flow  modelling  of  this  lower  deck  that  constitutes  the  generalized  global 
wall  layer  concept  of  this  research. 
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a)  Ue  conpleted  the  formulation  of  the  governing  Integral-type  equations 

that  characterize  the  non-reversed  flow  portion  of  Interactive  wall  layer 
model  In  the  form  of  a  set  of  nonlinear  partial  first  order  differential 
equations.  Including  the  three-dimensional  aspects  of  the  flow  (the 
cross  flow  velocity  component)  via  lateral  momentum  balance  :  see 

appendix  A. 

b)  A  significant  aspect  of  the  wall  layer  theory  is  the  further  splitting 
of  the  layer  into  an  upper  region  of  non-reversed  flow  and  an  underlying 
lower  region  of  slow,  reversed  flow  (see  Fig  3),  the  two  regions  being 

divided  by  a  locus  yu(x,z)  whose  unknown  shape  is  determined  from  the 

upper/lower  shear  stress  and  normal  velocity  matching  across  it. 
Incorporation  of  this  bi-regional  aspect  of  the  wall  layer  provides  the 
proper  Interactive  viscous/turbulent  physics  that  underlies  regions  of  3-D 
SBLI-lnduced  separation.  We  have  also  completed  the  formulation  of  the 
3-D  equations  that  govern  this  aspect  of  our  interactive  wall  layer  model 
(Appendix  A)  along  with  a  paper  describing  their  validation  and  application 
in  the  case  of  laminar  2-D  flow  (Appendix  B) 

c)  An  Important  feature  is  to  capture  the  shock  structure  and  nonlinear 
rotational  inviscld  flow  behaviour  of  the  overlying  region  by  means  of  a 
suitable  3-D  shock/boundary  layer  applications  which  include  separation. 
After  many  discussions,  we  established  a  collaboration  with  Dr  Peter 
Gnoffo  at  NASA  Langley  to  use  his  "Laura"  codef  This  code  was  made 
operational  at  ISU  including  its  adaptation  to  accept  the  matching 
conditions  that  allow  the  aforementioned  wall  layer  equations  to  be 
Inserted  as  an  "interactive  overlay"  including  a  User  Manual  we 
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. . . . 

have  written  for  this  ISU  version  with  the  interactive  underlay 
modification  aspects  (Appendix  C). 

Regarding  this  part  of  the  work,  several  unique  aspects  of  our  wall  layer 

underlay  approach  are  noteworthy.  First,  by  virtue  of  our  previous 

7-/6 

experience  with  3-D  interactive  flows  we  can  provide 

particularly  effective  initial  conditions  that  are  analytically  derived 
from  steady  triple  deck  theory,  rather  than  just  arbitrary  numbers, 
thereby  significantly  accelerating  time-convergence  of  the  solution. 
Second,  we  would  have  been  able  to  independently  check  out  the  3-D  Euler  code 
solution  convergence  in  the  regime  of  weaker  shock  strengths  at  lower 
supersonic  external  flow  Mach  numbers  by  means  of  the  analytical  3-D  small 
disturbance  middle-deck  theory  developed  previously  by  the  author. 
Third,  (see  Fig.  4)  we  are  able  to  avoid  the  near  wall  difficulties 
experienced  by  Euler  codes  within  nonuniform  flow  regions  such  as  a 
boundary  layer  profile;  using  a  wall  layer  underlay,  such  codes  need  not 
be  solved  all  the  way  to  the  floor,  but  only  down  to  some  cutoff  height 
h(x,z)>0  where  u(x,h,z)  and  w(x,h,z)  are  non-zero.  We  further  avoid  a 
spurious  normal  mass,  momentum  and  energy  fluxes  into  the  flow  at  this 
height  by  not  choosing  fl  arbitrarily  but  rather  equal  to  the  interactive 
displacement  thickness  y^(x,z)  of  the  interactive  wall  layer  that  actually 
underlies  the  middle  deck.  This  manner  of  incorporating  the 
effect  of  the  wall  region  shear  disturbance  physics  not  only  provides 
computational  advantages  but  is  unique  to  our  approach. 

d)  Having  completed  the  detailed  theoretical  formulation  of  the  wall  layer 
model  and  its  combination  with  a  code  for  the  overlying  middle/outer  deck 
rotational-flow  regions,  the  resulting  united  programme  had  begun  to  be 
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applied  to  the  analysis  of  several  specific  interacting  flow  problems  for 
which  extensive  experimental  data  and  OT)  Navier-Stokes  numerical 
solutions  results  are  available  for  validating  comparisons.  As  a  first 
step,  we  completed  the  calculation  of  a  set  of  interactive  solutions  for 
adiabatic  2-d  laminar  flow  past  a  compression  corner  that  include 
substantial  separation  (see  Appendix  D);  the  results  for  local  skin 
friction  and  pressure,  which  provide  a  significant  test  of  the  coupling  of 
our  wall  layer  model  and  its  ability  to  handle  reversed  flow,  are  in 
excellent  agreement  with  the  Navier-Stokes  solutions  of  Rizzetta.  This 
work,  the  last  to  be  completed  before  Grant  expiration,  is  being  written  up 
as  a  technical  paper. 
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Appendix  A 

Summary  of  the  Governing  Equations  characterizing  the  Interactive  Wall 

Layer  Model  as  a  CFD  Underlay. 

To  bring  out  the  physical  essentials,  we  further  adopt  the  following 
assumptions.  (1)  The  incoming  undisturbed  turbulent  boundary  layer  is 
characterized  by  a  steady  compressible  two-dimensional  Law  of  the  Wall/Law 
of  the  Wake  structure.  (2)  The  turbulent  mean  flow  within  the  interactive 
wall  sublayer  obeys  Morkovins  Hypothesis  as  regards  compressibility 
effects  and  is  modeled  in  the  first  approximation  by  Law  of  the  Wall 
mixing  length  concept  (3).  The  laminar  and  turbulent  Prandtl  numbers  are 
both  unity.  (4)  The  normal  pressure  gradient  (dp/dy)  and  non-boundary 
layer  type  stress  terms  are  negligible  across  the  Interactive  wall  layer. 

Non  -Reversed  Flow  Region 

Under  the  aforementioned  assumptions,  the  3-D  interactive  flow  in  the 
non-separated  portion  of  the  wall  layer  is  governed  by  the  following  set 
of  equations: 
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(10) 

(11) 


(12a) 


(12b) 


(12c) 


i=Al  y[l-EXP(-y/A)] 


(12d) 


(12e) 


Relations  (12)  represent  Rotta’s  non~isotropic  Law  of  the  Wall  turbulent 
eddy  viscosity  model^*  which  includes  both  pressure  gradient  effects  and 
a  non-isotropy  factor  T^Cwhere  Tj^  =  I  for  the  isotropic  case).  For  non- 

adiabatic  wall  conditions,  especially  at  hypersonic  external  Mach  numbers, 
we  retain  (unlike  classical  asymptotic  deck  theory)  the  variable  density 
effect  across  the  inner  deck.  The  solutions  of  Eqs.  6-12  must  satisfy 
boundary  conditions  with  respect  to  y  that  the  disturbance  shear  stresses 
vanish  beyond  some  inner  sublayer  thickness  5gj^(x,z)  and  that  on  y  =  0  we 

have  the  impermeable  floor  no-slip  boundary  conditions  V  (x,z)  =  U  (x,z)  = 
W(x,z)  =  0.  In  addition,  we  impose  the  upstream  condition  that  the 
flow  pass  over  to  the  given  compressible  turbulent  Law  of  the  Wall/Wake 
boundary  layer  profile;  on  the  basis  of  the  Reference  Temperature  version 
of  Waltz’s  formulation,  for  example,  these  profiles  in  the  general 
non-adlabatlc  wall  case  are  described  analytically  by 


To(y)  -  T^+  (Tj 


Where  t)  =  y/5o,  R  =  4.1  Re5o/((l  +  n)  (Tw/Te)'^*^]  with  u  =  .76  and  n  is 
Coles’  wake  function  linked  to  Cfo  and  Redo  by  the  relation 
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(15) 


27C  +  .215  +  ln(l+R)  =  .41/ 

Equation  (13)  for  tj  >.10  yields  a  Law  of  the  Wake  behaviour  with  UoAJe  -»  1 
and  dUo/dy  ->  0  as  t)  =  1,  while  for  very  small  v  it  assumes  a  Law  of  the 
Wall-type  behavior  consisting  of  a  logarithmic  term  that  is  exponentially 
damped  out  into  the  linear  sublayer  profile  u/u  =  Rt)  as  i)  o. 


Reversed  Flow  Region 


Turning  to  a  consideration  of  any  reversed  disturbance  flow  that  might 
occur  in  the  lower  region  of  the  wall  layer,  we  propose  to  treat  it  on  the 
basis  of  a  negligible  inertial  approximation  that  is  appropriate  to  such 
slow  wall-dominated  flows.  Thus,  we  must  solve  the  following  equations  by 
a  y-integration  of  Eqs.  S-11  with  the  convective  acceleration  terms 
dropped: 


3u 

n— +pvp„=T  + 
By 


9Pw 


(16) 


Bv/ 

By 


(17) 


(18) 


Where  we  have  used  the  fact  that  vanishes  at  the  wall.  Equations  16-18 
are  to  be  solved  up  to  the  vertical  height  yu(x,z)  of  the  Un  =  0  locus 
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(unknown  priori)  where  Un  «  U  sin  A  -  U  cos  A  and  A  is  the  angle  of  the 
local  separation  line  as  shown  in  Fig  5. 


Outline  of  the  Wall  Solution  Method 

Addressing  first  the  reversed  flow  region-aspect  of  this  wall  layer  for 

0  s  y  s  yu  (x,z),  the  velocity  and  temperature  (and  density)  field 

solutions  are  to  be  obtained  by  further  integrations  wrt  y  of  Eqs.  16-19 

along  with  equations  12a-e.  When  the  resulting  expression  for  Un  is 

equated  to  zero,  we  additionally  get  a  pair  of  expressions  governing  the 

unknown  locus  Yu(x,z)  of  the  form 

^yu  ^Pw  ,  ,  , 

— +  J  (19A) 

dx  ax 

+  •  ^z,  'z^yu)  (19B) 

oz  oz 

where  the  various  K’ s  here  are  detailed  functionals  resulting  from  the 
integration.  An  additional  integration  of  the  continuity  Equation  then 
yields  an  expression  for  the  transpiration  mass  flux  across 


U  of  the  form 


y,(*.b) 


PvCx.Yu.z)  = 


3(pu)  ,  3(pw)  L 
- 1 - ay 

dx  dz 


2  ^P 

yu  ~r  ^"^x j^z  'yu) 

dx 


yu  iz3  ('^xj'^z.'yu) 

dz 

Where  the  J’s  are  further  integration-functionals.  Along  with  us,  we  note 
that  a  further  y-lntegratlon  of  Eqs.  16-18  yield  relationships  between  the 
wall  shear  components  and  the  shears  along  the  in  terms  of  8p/dx  and 
dp/dz  of  the  form 
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layer  for  purposes  of  a  CFD  underlay,  we  obtain  a  solution  of  the 
equations  by  an  integral  method  approach:  this  not  only  enables 
an  efficient  matching  between  the  upper  and  lower  wall  layer  regions 
across  Yu  but  also  formulates  the  wall  layer  model  in  the  desired  form  of 
a  set  of  first  order  partial  differential  equations  governing  the  3-D 
disturbance  wall  shear,  heat  transfer  and  displacement  thickness  Yw(x,z). 
Indeed,  some  unpublished  earlier  work  by  R  T  Davis  reveals  that  it  yields 
relationships  of  the  form 


(26) 

(27) 

(28) 


where  the  various  Fx,  Fz  and  Fo  here  are  functionals  of  Tx(yu),  pH(x,y) 
and  yu(x,z)  given  by  the  details  of  the  momentum  and  energy  equation 
integral  analyses  with  suitable  profiles.  When  combined  with  Eqs.  19-25 
that  link  Tx(yu),  Tz(yu)  and  q(yu)  along  y  «  yu  to  their  local 
counterparts  along  the  wall  y  =  0,  Eqs.  26-28  provide  in 
principle  a  trio  of  relationships  governing  the  local  interactive  skin 


(29) 

(30) 


(31) 


with  the  aforementioned  analysis  of  the  entire  wall  layer  completed,  we 
close  the  problem  by  obtaining  the  overall  vertical  velocity  field  from  an 
integration  of  continuity  Eq  6  across  0  s  y  s  yu  and  then  beyond  the 
viscous  thickness  of  the  wall  layer  6sl  where  the  disturbauice  shear 
stresses  become  vanishingly  small;  this  procedure  yields 


When  the  right  side  of  this  Equation  is  equated  to  zero,  the  result 
provides  a  determination  of  the  y(x.z)  locus  along  which  the  "far  field"  v 
emerging  from  the  wall  layer  vanishes,  ie,  the  effective  inviscid  wall 
position  (or  displacement  thickness)  yw(x,z)  (see  again  fig  4.),  The 
resulting  expression  for  yw  that  emerges  is  of  the  general  form. 


ax  ”  az  (33) 

Eq.  33  completes  the  formulation  by  providing  the  desired  coupling 
mechanism  to  the  overlaying  numerical  Euler  solution  of  the  middle  deck. 
Moreover,  the  detailed  evaluation  of  Eq.  32.  at  the  outer  edge  of  the 
proposed  generalized  wall  layer  also  permits  a  careful  examination  of  the 
important  property  of  the  upwash  (vertical  velocity)  behaviour  right  above 
3-D  skin  friction  line  confluence.  This  is  very  helpful  in  further 
illuminating  the  actual  physical  behavior  of  the  flow  as  it  ostensibly 
leaves  the  surface  ("separates")  along  such  a  confluence. 
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1  INTRODUCTION 

'Jin'  |ir()^i';iii)  I.AHUA  u.is  \vnii>'ii  .il  NASA  l,:ui};li‘y  l>y  I’l'IiT  (inollu.  Ait> 
(ini's)  ions  iiol  .iiisni'H'i!  l>y  iliis  i  :in  Ix'  iliri'Cli'il  lo  Itiin  nl  liis  l>iiiail 

aililii'ss  III i(l< 1 1; .iIk  I' 111  I  l.i|(  nasa.s;o\  ‘  or  liy  ti'li'jilioiii*  at  8()-l-‘l-J80.  lli*  is 

vi'iy  ln'l|>rnl  ami  i.-'iy  (Mli.-ni.  I’l-l.-i  ami  I  iliil  inosi  of onr  inli'rarlion  llir')nji;|i  l-'.- 
iiiail.  I  sent  liiin  i|m-slioiis.  ami  In' . soul  ino  answers  ami  jirograiii  iiiotliliral.iuiis. 

'I'll*'  roili-  ilsi-lf  is  ^I'ly  larg*',  ami  coiisisJs  of  aUonI  91)  Sr»|iarnl€*  I'urlran 

lili's,  I'l'li'i  inaili'il  ili.'s,'  |i|(si  lo  my  I’lojct'l.  Viiicenl  affoiint,  ami  we  wi-til 
to  work  Iryinjj  lo  j;'-!  a  vvoikinj*  versicni  lii'rc  at  ISl’.  (Ielliii|'  tin*  prograiti  lo 
compile  took  several  u.-i  ks  of  inleractioii  since  it  liad  originally  Inx'ii  conipilerl 
oil  a  dilferenl  maclilm;.  ami  a  few  of  tlie  Ofttions  liatl  lx.vn  slrippeil  (i.c.  reurtiiig 
flow  ami  iiuill ilasking  capaliilitics).  Finally,  however,  we  tiiti  gel  the  program 
to  work. 

Although  I  lie  progiam  is  capahle  of  modeling  3-D,  viscous,  sn  person  ie/hy  person  ic 
How  in  tliermocln-mical  iioneiiiiilihrinm,  we  will  be  primnrilly  interested  in  2-D, 
inviscid  sui>ersoiiic,  noiiieacting  How  with  vorlicity. 

2  INPUT 

2.1  wingbdy.f 

'riie  geomolry  of  the  liody  is  containe*!  in  the  file  “winghdy.P.  I'hc  exam|>le 
body  that  is  stored  there  now  is  a  2-U  cylindrially  blunted  15  degree  wialge, 
which  was  chasen  .so  that  the  shock  will  be  curved,  and  the  How  will  have 
vorticity. 

2.2  Input  File 

The  the  input  file  is  used  to  cnstomir,e  the  flow  and  run  parameters.  An  e.\artiip|e 
iiipiil  file  is  shown  in  riguie  1.  'I'he  variable  iiaiiics  and  their  significance  is  given 

• 

# 

_ _  • 

5.3  Inger.q 

Dr.  fttger  is  doing  work  that  requires  cerlaiii  informal ir>n  at  or  near  the  surface 
of  the  body.  Tiierfore,  Peter  has  modified  l,AUHy\  lo  produce  a  formatted 
output  file  containing  tliis  information  called  “ingern". 

The  first  row  of  this  file  contains  a  single  mmilx'r  ''I  '.  'I’his  di'noti's  that 
the  block  of  information  that  follows  is  at  the  hotly  surface.  Several  lines  down 
from  this  single  ‘'1"  is  another  row  tliat  contains  a  siiiKle  number  “2".  'I  his  “V” 
denotes  that  the  block  of  information  that  follows  is  for  the  cells  that  are  a  step 
away  from  the  body.  After  the  "2"  block,  there  is  a  “.I”  block  that  is  one  more 
step  away  from  the.  body.  Uy  selecting  the  proper  block,  the  us«*r  can  use  the 
fiow  properties  at  the  body  surface  or  somwhere  above  (he  hotly  surface. 

Each  block  of  information  contains  several  rows  of  infurmalion.  Each  row 
has  the  following  eleven  elements: 

-  ~  ^  p  u.  i!u  ie.  iiU.iX- 
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ABSTRACT 

Separation  Bubble  Modeling  in  Compression 
Ramp-Generated  Viscous-Inviscid  Interactions 

G.  R.  Inger  and  B.  Nath* 

Iowa  State  University 

iNTRODUCnON 

Correct  treatment  of  the  near  wall  viscous  flow  physics,  and  especially  the  separation  bubble 
region,  is  an  important  aspect  of  the  calculation  of  the  interacting  flow  field  generated  around  a 
deflected  control  surface  on  a  high  speed  aerodynamic  body.  Such  near  wall  features  are  often 
difficult  and  expensive  to  resolve  with  a  CFD/Navier-Stokes  code  without  fine  grained  sub-mesh 
refinement  schemes.  It  would  therefore  be  of  value  in  practical  design-type  codes  if  an  analytical 
underlay  to  them  could  be  supplied  that  resolves  the  near  wall  viscous  flow  physics  including  the 
reversed  flow  ponion  of  the  separation  bubble.  This  paper  presents  such  an  underlay  theory  for 
the  simplified  trial  problem  of  2-D  steady  laminar  supersonic  flow  past  a  compression  comer. 
OUTLINE  OF  THE  THEORETICAL  APPROACH 

Consider  the  vicinity  of  a  high  Reynolds  number  compression  ramp  flow  as  schematically 
illustrated  in  Fig.  1:  here,  as  is  well-known,^  the  local  strong  viscous-inviscid  interaction 
generates  a  triple  deck  structure  outside  and  inside  the  incoming  boundary  layer  profile,  consisting 
of  an  outer  region  of  disturbed  irrotational  inviscid  flow,  a  middle  deck  of  rotational  perturbation 
flow,  and  an  inner  deck  of  non-reversed  viscous  disturbance  flow.  This  triple  deck  in  turn  can  be 
imagined  to  lie  on  top  of  the  portion  of  the  separation  bubble  that  contains  a  revmed  (viscous) 
flow. 

Our  treatnnent  of  this  situation  consists  of  two  parts.  In  the  overlying  triple  deck  of  non- 
reversed  flow,  we  have  developed  a  unique  defect-type  of  integral  method  which  describes  the 
overall  local  features  of  the  deck  in  terms  of  a  set  of  ordinary  differential  equations.  As  shown  in 
Thus  if  pvvr,  5,  Ui  denote  the  local  wall  pressure,  thickness  and  edge  slip  velocity. 

Professor  and  graduate  student,  respectively,  Dept,  of  Aerospace  Engineering  and  Engineering 
Mechanics 


respectively  of  the  deck,  these  equations  for  supersonic  outer  deck  flow  are  as  follows  in  terms  of 
the  well-known  non-dimensional  triple  deck  variables  (Ref.  1): 

where  a  is  the  ramp  angle,  Yu(x)  is  the  underlying  separation  bubble  locus  upon  which  this  deck 
rests  (see  Figjjand  Vu  is  the  corresponding  transpiration  velocity  across  Yu(x)  into  the  bottom  of 
the  deck  from  the  reversed  How  region.  It  is  seen  that  the  solution  of  Eqs.  (l)-(3)  is  necessarily 
coupled  to  the  behavior  (unkno'.'.  n  apriori)  of  the  reversed  separation  bubble  properties. 

The  second  aspect  of  our  theory  is  thus  to  provide  a  comparable  analytical  treatment  of  the 
reversed  flow  zone  whereby  its  Yu(x)  and  Vu(x)  arc  appropriately  linked  to  the  pTq>erties  of  the 
overlying  triple  deck.  This  has  been  done  using  a  method  of  successive  inertial  iq)proximations;^ 
accord!  ng  to  the  detailed  analysis  to  the  full  paper,  the  results  are  the  following  : 

5^=0  if  2Vl/S<l  1 

V.- 

The  simultaneous  solution  of  Eqs.  ( l)-(5)  for  the  five  unknowns  p,  Ui,  5,  Vu  and  Yu  comprises 
the  matched  treatment  of  the  entire  interaction  zone  including  the  reversed  flow  aspect  of  the 
separation  bubble  that  occurs  for  large  enough  compression  comer  angles. 

VALIDATION  OF  THE  THEORY 

The  reversed  flow  ponion  of  the  foregoing  theory  has  been  validated  by  comparison  with  the 
now  classical  exact  Navier-Stokes  numerical  solutions  by  Briley  for  a  typical  separation  bubble 
flow.3  As  shown  in  Fig.  2,  the  present  theory  gives  very  accurate  results  fw  the  Yu  kxnis  in  such 


a  flow,  and  hence  should  provide  a  faithful  representation  of  die  reversed  flow  aspect  of  the  above- 
described  interaction  thetHy. 

The  foregoing  triple  deck  theory  of  the  reversed  flow  region  also  has  been  shown  to  be  in 
good  agreement  with  exact  numerical  solutions  :  an  example  is  illustrated  in  Fig.  3, 
(n  the  case  Si  =  2.50. 

TYPICAL  RESULTS 

In  Figure  4  are  illustrated  a  set  of  representative  results  obtained  by  a  Runga-Kutta  method 
solution  of  Eqs.  (l)-(5)  for  a  separated  flow  case  with  a  =  2.2.  The  significant  local  alterations  of 
the  interactive  pressure  and  skin  friction  by  the  comer-smoothing  displacement  effect  of  the  bubble 
are  closely  seen.  Currently  this  work  is  being  extended  to  address:  (a)  hypersonic 
external  flow,  (b)  non-adiabatic  wall  conditions  with  heat  transfer  across  the  separation  bubble,  (c) 
three  dimensional  flows  and  (d)  turbulent  flow  including  nirbulence  modeling  within  the  reversed 
flow  region. 
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Fig.  1  MuItiplc-Dcckcd  Configuration  of  Compression  Ccnner  Interaction  Regitxi 


Fig.  3  Comparison  of  Present  Interaction  Solution  with 
Separation  with  Exact  Solutions  of  Rizzetta 
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Figures  for  Report 


Reversed  Flov/  Region 

Fig.  3  Detail  of  Reversed  Flow  Structure  at  Bottom  of  Wall  Layer 
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Fig.  4  Middle  Deck  Matching  to  the  Wall  Layer 
by  the  Displacement  Thickness  Concept 


